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A  simple  circuit  model  is  developed  to  describe  supercapacitor  behavior,  which  uses  two 
resistor-capacitor  branches  with  different  time  constants  to  characterize  the  charging  and  redistribution 
processes,  and  a  variable  leakage  resistance  to  characterize  the  self-discharge  process.  The  parameter  val¬ 
ues  of  a  supercapacitor  can  be  determined  by  a  charging-redistribution  experiment  and  a  self-discharge 
experiment.  The  modeling  and  characterization  procedures  are  illustrated  using  a  22F  supercapacitor.  The 
accuracy  of  the  model  is  compared  with  that  of  other  models  often  used  in  power  electronics  applications. 
The  results  show  that  the  proposed  model  has  better  accuracy  in  characterizing  the  self-discharge  process 
while  maintaining  similar  performance  as  other  models  during  charging  and  redistribution  processes. 
Additionally,  the  proposed  model  is  evaluated  in  a  simplified  energy  storage  system  for  self-powered 
wireless  sensors.  The  model  performance  is  compared  with  that  of  a  commonly  used  energy  recursive 
equation  (ERE)  model.  The  results  demonstrate  that  the  proposed  model  can  predict  the  evolution  profile 
of  voltage  across  the  supercapacitor  more  accurately  than  the  ERE  model,  and  therefore  provides  a  better 
alternative  for  supporting  research  on  storage  system  design  and  power  management  for  wireless  sensor 
networks. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Harvesting  energy  from  ambient  environment  to  power  wire¬ 
less  sensor  networks  (WSNs)  has  been  investigated  to  extend 
system  lifetime  [1-8].  Rechargeable  batteries  such  as  NiMH  [1] 
and  Li-ion  [2]  were  first  selected  to  serve  as  primary  energy  stor¬ 
age  devices.  While  rechargeable  batteries  have  high  capacity  and 
low  leakage,  the  cycle  life  of  rechargeable  batteries  limits  the 
lifetime  of  wireless  sensors  [3,5,6].  The  cycle  life  of  a  recharge¬ 
able  battery  is  defined  as  the  number  of  charge-discharge  cycles 
before  its  capacity  falls  below  80%  of  its  initial  rated  capacity. 
The  aging  process  during  the  charge-discharge  cycles  results  in 
a  gradual  reduction  in  capacity  and  an  increase  in  internal  resis¬ 
tance  over  time  [9].  By  the  end  of  the  cycle  life,  the  capacity  of 
a  rechargeable  battery  is  reduced  by  20%  and  the  useful  energy 
drops  to  50%  because  the  higher  internal  resistance  causes  prema¬ 
ture  end  of  life  [3].  Due  to  limited  cycle  life,  typically  ranges  from 
100  s  to  1000  s  [10],  a  wireless  sensor  node  will  require  battery 
replacement  after  1-2  years  [3,5,6].  Compared  with  rechargeable 
batteries,  supercapacitors  have  a  much  longer  cycle  life  and  higher 
charge-discharge  efficiency  in  addition  to  fast  charge-discharge 
characteristic  [3,5-7,11,12].  Some  authors  have  proposed  to  use 
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supercapacitors  alone  to  store  harvested  energy  [3,12]  or  use  super¬ 
capacitors  in  combination  with  rechargeable  batteries  [5-7,11]  to 
achieve  “perpetual  lifetimes”  for  wireless  sensor  networks. 

Due  to  wide  range  of  potential  applications  of  self-powered 
WSNs,  there  is  considerable  research  interest  in  developing 
efficient  power  storage  systems  [11,12],  power  management  algo¬ 
rithms  [13]  and  communication  protocols  [14].  A  supercapacitor 
model  is  an  important  tool  for  evaluating  these  researches  using 
analytical  methods  or  simulation  before  being  demonstrated  in 
practical  deployments.  Currently  adopted  supercapacitor  models 
in  WSN  research  are  developed  from  the  supercapacitor  leakage 
power  profiles  [5,11,12,14,15],  which  are  called  the  energy  recur¬ 
sive  equation  (ERE)  models  in  this  paper.  The  ERE  model  assumes 
the  capacitance  of  a  supercapacitor  is  constant,  and  uses  a  long¬ 
term  leakage  power  profile  to  determine  leakage  at  any  time  [12]. 
However,  the  physics  of  the  supercapacitor  suggests  that  its  capac¬ 
itance  depends  on  the  terminal  voltage  across  the  device  and  the 
supercapacitor  usually  experiences  internal  charge  redistribution 
during  and  after  charge  and  discharge  cycles  [16].  Without  captur¬ 
ing  these  properties,  the  terminal  voltage  estimated  using  the  ERE 
model  may  have  significant  deviation  from  the  actual  value.  The 
voltage  deviation  can  cause  improper  decision  making  in  a  power 
management  system  for  self-powered  wireless  sensors  that  harvest 
ambient  energy  and  experience  frequent  charge-discharge  cycles. 
The  improper  decision  may  even  lead  to  failure  of  a  wireless  sensor 
network. 
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A  simple  equivalent  circuit  model  for  supercapacitors  that  can 
accurately  model  charging,  redistribution  and  self-discharge  pro¬ 
cesses  is  presented  in  this  paper.  A  supercapacitor  with  rated 
capacitance  of  tens  Farads,  typical  capacity  for  WSN  applications,  is 
used  as  an  example  to  illustrate  the  modeling  and  characterization 
process.  The  performance  of  the  developed  model  is  compared  with 
that  of  models  used  in  power  electronics  applications  that  often 
target  supercapacitors  with  rated  capacitances  of  kilo-Farads.  In 
addition,  the  developed  model  is  compared  with  the  ERE  model  in 
analyzing  a  simplified  energy  storage  system  for  wireless  sensors. 
The  results  demonstrate  that  the  new  model  provides  a  more  accu¬ 
rate  estimation  of  terminal  voltages  across  the  supercapacitor  than 
the  currently  used  ERE  model  and  hence  a  better  alternative  for  sup¬ 
porting  research  on  storage  system  design  and  power  management 
in  WSNs. 


Time  (s) 

Fig.  2.  A  22F  supercapacitor  charged  with  a  constant  97  mA  current. 


tion  during  charge  and  discharge  cycles  [16,31  ],  The  second  branch 
(R2  and  C2)  is  the  delayed  branch,  which  represents  the  medium 
and  long  term  charge  redistribution  [16,31  ].  The  voltage  relaxation 
is  the  result  of  charge  redistribution  [32].  The  equivalent  parallel 
resistance  (EPR)  represents  the  effect  of  self-discharge  [16,31  ]. 

2.1.  R0 


2.  Modeling  and  characterization  of  supercapacitors 

The  electrochemical  impedance  spectroscopy  (EIS)  is  a  general 
approach  to  measure  the  complex  impedance  of  energy  storage 
devices  such  as  supercapacitors  and  batteries  [17].  The  nature  of 
impedances  in  various  frequency  ranges  can  be  determined  by 
analyzing  the  frequency  dependencies  of  the  real  part  and  the  imag¬ 
inary  part  [18].  Various  equivalent  circuit  models  [  1 9,20]  have  been 
developed  using  the  porous  electrode  theory  [21-26]  to  interpret 
the  impedance  spectrum  of  a  supercapacitor.  With  the  assumption 
of  homogeneous  electrode  pore  size,  a  general  porous  impedance 
model  consists  of  three  impedances  linking  to  electrode,  electrolyte 
and  electrode/electrolyte  interface  [23-25].  The  general  model 
can  be  modified  if  interface  roughness  [27]  or  random  pore  size 
[28]  is  considered.  The  impedance  spectrum  can  be  modeled  by 
N  interleaved  RC  circuits  which  would  require  the  calculation  of 
2N  parameters  [17].  However,  it  is  usually  very  hard  to  determine 
more  than  five  or  six  independent  parameters  efficiently  due  to  the 
strong  influences  between  them  [17]. 

Alternatively,  supercapacitors  can  be  characterized  in  time 
domain  by  conducting  various  experiments  such  as  constant  power 
tests  and  constant  current  tests  [18].  In  power  electronics  applica¬ 
tions,  this  approach  is  often  used  to  develop  an  equivalent  circuit  to 
model  the  terminal  behavior  of  supercapacitors  [  1 6,29].  As  superca¬ 
pacitors  are  usually  constructed  with  two  porous  carbon  electrodes 
impregnated  with  electrolyte  and  separated  by  a  porous  insulat¬ 
ing  membrane,  the  interface  electrochemistry  suggests  a  complex 
network  of  non-linear  capacitors  connected  between  them  by  resis¬ 
tances  can  describe  the  behavior  of  supercapacitors  [29].  In  the 
network,  resistors  represent  the  resistivity  of  carbon  particles  and 
capacitors  describe  the  capacitance  between  carbon  and  electrolyte 
[16].  In  practical  applications,  simplified  versions  are  required  to 
enable  model  implementation  [29].  One  often  used  version  is  the 
two-branch  model  [31  ]  shown  in  Fig.  1.  The  first  branch  is  the  main 
branch,  which  includes  R0  and  the  voltage  dependent  capacitance 
(C0  and  I<vV).  The  main  branch,  also  called  the  immediate  branch, 
dominates  the  immediate  behavior  of  a  supercapacitor  in  response 
to  a  charge  action  and  captures  the  linear  dependence  of  capac¬ 
itance  on  the  terminal  voltage  in  the  practical  voltage  range  of 
the  device  [16].  This  branch  determines  short-term  voltage  evolu¬ 


The  resistance  R0  can  be  determined  by  measuring  the  potential 
difference  A  V  between  the  two  terminals  at  the  beginning  while 
charging  the  supercapacitor  with  a  large  constant  charging  cur¬ 
rent  Jc  [16,31  ].  The  ratio  of  initial  short-term  potential  change  and 
/c  determines  the  R0  value.  As  the  R0  value  identified  using  this 
approach  is  very  close  to  the  ESR  (equivalent  series  resistance) 
value  provided  by  the  supercapacitor  datasheet,  it  is  reasonable 
to  use  ESR  value  for  R0. 


2.2.  C0  and  I<v 


Assuming  all  the  charge  is  stored  in  the  first  branch  during 
the  charging  process,  the  current-voltage  relationship  across  the 
capacitor  of  the  first  branch  is  [31  ]: 


dq  dv 
dv  dt 


=  (c0  +  i<w)Tt 


(1) 


In  case  of  a  constant  charging  current  Jc ,  the  time-voltage  rela¬ 
tionship  at  the  supercapacitor  terminals  can  be  derived  from  (1 )  as 
follows: 


t=m=<*v+l«v>  (2) 

Given  two  data  points  (ti,Vi)  and  (t2,V2)  in  the  time- voltage 
charging  curve  with  a  constant  charging  current  (Fig.  2),  the  fol¬ 
lowing  equations  must  hold: 


h 

< 

£2 


The  values  for  Co  and  Kv  can  be  solved  from  (3): 


(3) 


tl_  _  Vi  t2  ~  tl  V2 
Vi  v22-v,v2 


Ky  =  2 


Vlf2  -  C|V2 

v^vl-v}v2 


(4) 
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Fig.  3.  Self-discharge  of  a  fully  charged  22F  supercapacitor. 


C2  and  R2 

After  the  supercapacitor  reaches  its  rated  voltage,  the  charging 
current  is  removed.  The  charge  stored  in  the  first  branch  charges 
the  second  branch  during  the  redistribution  process  [3 1  ].  If  the  time 
constant  of  the  second  branch  r2  =  R2C2  is  fixed,  we  can  assume  the 
voltages  of  both  branches  are  same  at  time  3r2.  The  value  r2  =  240 
s  can  be  used  as  a  starting  point  [31].  If  it  does  not  provide  suf¬ 
ficient  accuracy  for  the  model,  more  experimental  tests  can  be 
performed  to  adjust  the  value.  Once  the  time  constant  of  the  sec¬ 
ond  branch  is  determined,  the  voltage  of  the  supercapacitor  at  time 
3r2,  V2f,  during  the  redistribution  process  can  be  measured.  At  this 
point,  assuming  all  charge  is  stored  in  the  supercapacitor,  C2  can  be 
determined  from  the  following  equation: 

Qtot  =  IcTc  =  C2 V2{  +  (c0  +  ^ V2f)  V2{  (5) 

where  Tc  is  the  charging  time.  Then  the  value  of  resistor  R2  is: 

R2  =  p-  (6) 


0 

CD 

B 

o 

> 


2.3. 


3.  Modeling  supercapacitor  self- discharge 

A  charged  supercapacitor  is  in  a  state  of  higher  Gibbs  energy 
than  in  its  discharged  state  [33].  Therefore  a  thermodynamic  “driv¬ 
ing  force”  results  in  spontaneous  decline  of  Gibbs  energy  [33].  This 
decline  manifested  as  decay  in  supercapacitor  voltage  is  the  self¬ 
discharge.  The  rate  of  self-discharge,  which  usually  diminishes  with 
time,  actually  determines  the  shelf-life  of  supercapacitors  [34]. 
Fig.  3  shows  the  self-discharge  of  a  22F  supercapacitor.  The  open  cir¬ 
cuit  voltage  of  the  supercapacitor  is  measured  for  8  h  after  it  is  fully 
charged  to  its  rated  voltage  2.5  V  with  a  constant  voltage  source.  It 
is  obvious  that  the  self-discharge  rate  decreases  with  time. 

The  open-circuit  self-discharge  of  supecapacitors  must  take 
place  through  coupled  anodic  and  cathodic  processes  to  pass  par¬ 
asitic  currents  at  one  or  both  individual  electrodes  since  there  is 
no  external  circuit  through  which  discharge  can  pass  [33].  The 
self-discharge  can  be  ascribed  to  three  mechanisms  [30,34,35]: 
activation-controlled  Faradaic  processes  if  the  supercapacitor  is 
overcharged,  diffusion-controlled  Faradaic  redox  reactions  if  the 
supercapacitor  is  not  overcharged  but  some  depolarizing  impurity 
is  present,  and  leakage  current  if  the  supercapacitor  has  internal 
ohmic  leakage  pathways.  Under  the  normal  operation  conditions, 
both  diffusion-controlled  Faradaic  redox  reactions  and  internal 
ohmic  leakage  need  to  be  modeled. 

The  two-branch  model,  in  which  the  self-discharge  is  repre¬ 
sented  with  a  constant  equivalent  parallel  resistance  (EPR),  only 


R1  R2 


considers  the  internal  ohmic  leakage.  To  model  the  nonlinear  self¬ 
discharge,  Diab  et  al.  [36]  proposed  a  three-branch  model  (Fig.  4) 
in  which  the  RrCv  branch  was  added  to  represent  the  self-discharge 
due  to  diffusion-controlled  Faradic  redox  reactions. 

The  leakage  current  due  to  diffusion-controlled  Faradic  redox 
reactions  is  proportional  to  the  concentration  gradient  of  the  dif¬ 
fusible  redox  species.  The  time-dependence  of  the  concentration 
gradient  at  a  particular  distance  from  a  plane  electrode  is  gener¬ 
ally  inversely  proportional  to  the  square-root  of  time,  which  leads 
to  a  decrease  of  self-discharge  rate  with  time  [34].  We  propose  to 
use  a  variable  leakage  resistance  (VLR)  to  model  the  time  depen¬ 
dency  of  self-discharge  (Fig.  5).  We  do  not  distinguish  different 
self-discharge  mechanisms.  Rather,  we  introduce  a  single  VLR  to 
take  into  account  the  effects  of  various  self-discharge  mechanisms 
under  the  normal  operation  conditions. 

With  a  constant  leakage  resistance  EPR  (Fig.  1),  the  self¬ 
discharge  can  be  described  by  the  following  exponential  function 
[36]: 


f(t)  =  Vo  exp  (-T)  (7) 

where  V0  is  the  initial  voltage  and  X\e  is  the  time  constant  of  leakage 
current. 

Once  the  time  constant  of  leakage  current  is  determined  exper¬ 
imentally,  the  leakage  resistance  EPR  can  be  calculated  based  on 
the  following  equation  [36]: 


EPR  =  5^ 

M 


(8) 


with  a  variable  leakage  resistance,  VLR,  instead  of  using  one 
exponential  function  to  model  the  self-discharge  characteristic,  we 
use  several  distinct  exponential  functions.  These  exponential  func¬ 
tions  have  different  time  constants  during  different  self-discharge 
periods,  which  result  in  different  leakage  resistances.  By  trial  and 
error,  we  found  that  five  distinct  exponential  functions  (Eq.  (9)) 
corresponding  to  five  periods  can  fit  the  self-discharge  measure¬ 
ment  of  a  supercapacitor  up  to  28,800  s  (8h)  very  well  (Fig.  6(a)). 
Fig.  6(b)  shows  a  zoomed-in  result  for  the  self-discharge  during  the 
first  2  h.  For  the  following  sections,  the  self-discharge  characteristic 
will  be  restricted  to  the  first  2  h  because  we  are  more  interested  in 
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Fig.  6.  Comparisons  between  self-discharge  measurements  and  fitted  multiple 
exponential  functions:  (a)  8  h  and  (b)  2  h. 


the  short  and  medium  term  behavior  of  supercapacitors. 


v(t)  = 


2.504  exp(-1.181  x 
2.488  exp(-7.807  x 
2.462  exp(-5.301  x 
2.425  exp(-3.279  x 
2.382  exp( -2.042  x 


10“5t), 

t  =  (0,1500  s) 

10“6t), 

t  =  (1500,  4000  s) 

10“6t), 

t  =  (4000,  7200  s) 

10“6t), 

t  =  (7200,  14400  s) 

10“6t), 

t  =  (14400,  28800  s) 

(9) 


4.  Evaluation  of  the  VLR  model 

The  VLR  model  for  a  22F  supercapacitor  from  Cooper  Buss- 
mann  [37]  was  developed  following  the  characterization  process 
described  in  Sections  2  and  3.  The  charging,  redistribution  and 
self-discharge  processes  were  performed  for  identifying  model 
parameter  values.  During  the  charging  and  redistribution  pro¬ 
cesses,  the  initially  completely  discharged  supercapacitor  was 
charged  for  610  s  using  a  constant  current  source  (97  mA).  Then  the 
current  source  was  turned  off  and  the  decline  of  voltage  across  the 
supercapacitor  due  to  redistribution  was  tallied  for  the  next  890  s. 
During  the  self-discharge  process,  we  assumed  that  the  superca¬ 
pacitor  was  fully  charged  after  being  charged  for  1  h  with  a  2.5  V 
constant  voltage  source.  The  self-discharge  was  recorded  for  2  h 
after  removing  the  voltage  source. 

To  reduce  random  measurement  errors,  each  experiment  was 
repeated  three  times  to  obtain  an  average.  Fig.  7  shows  the  experi¬ 
mental  data  and  the  corresponding  average  values.  The  root  mean 


Time  (s) 

Fig.  7.  Comparisons  between  three  experiments:  (a)  charging  and  redistribution 
processes  and  (b)  self-discharge  process. 


square  error  (RMSE)  between  the  measured  experimental  data  and 
the  average  is  below  0.05  V  for  the  charging  and  redistribution  pro¬ 
cesses  (Fig.  7(a)),  and  below  0.005  V  for  the  self-discharge  process 
(Fig.  7(b)).  The  RMSEs  are  relatively  small  compared  to  the  voltage 
ranges  involved  in  the  experiments.  The  model  parameter  values 
of  the  22F  supercapacitor  (Table  1 )  were  determined  based  on  the 
average  data.  The  simulation  results  of  the  model  match  the  aver¬ 
age  of  measurements  very  well  for  charging,  redistribution  and 
self-discharge  processes  as  shown  in  Fig.  8. 

Performance  of  the  proposed  VLR  model  is  compared  with  that 
of  two-branch  and  three-branch  models.  The  constant  leakage 
resistance  EPR  in  a  two-branch  model  can  be  determined  using  the 
ratio  of  the  nominal  voltage  and  the  leakage  current  reported  in  the 
datasheet  [31].  For  the  22F  test  sample,  the  leakage  current  is  not 
specified  in  the  manufacturer  datasheet  [37].  To  estimate  the  EPR, 
we  used  an  exponential  function  to  fit  the  self-discharge  curve  [36] 
and  found  a  leakage  resistance  value  for  EPR  (4350  £2)  based  on  Eqs. 
(7)  and  (8).  A  larger  value  (9000  Q)  and  a  smaller  value  (2000  Q) 
were  also  used  for  comparison. 

The  comparison  of  different  models  is  shown  in  Fig.  9.  All  mod¬ 
els  have  almost  the  same  performance  during  the  charging  process 
(Fig.  9(a)),  and  can  accurately  predict  the  voltage  values  across 


Table  1 

VLR  model  parameters  of  a  22F  supercapacitor  from  Cooper  Bussmann. 


Model  parameter 

Co 

I<v 

Ro  [37] 

c2 

r2 

Variable  leakage  resistance  (VLR) 

0-1500 s  1 500-4000 s 

4000-7200 s 

Value 

16.2F 

4.8F/V 

0.04  Q, 

2.3F 

210Q 

2956  Q  4398  Q 

6523  Q, 
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Fig.  8.  Comparisons  between  the  experiment  and  the  simulation  results  of  VLR 
model:  (a)  charging  and  redistribution  processes  and  (b)  self-discharge  process. 


the  supercapactior  during  the  linear  charging  region.  At  the  end 
of  the  charging  process,  the  voltage  value  is  slightly  underesti¬ 
mated  by  all  models.  During  the  redistribution  process  (Fig.  9(a)), 
three-branch  model  has  the  best  performance,  but  other  models 
also  have  fairly  accurate  prediction  except  the  two-branch  model 
with  EPR  of  9000  £2.  During  the  self-discharge  process  (Fig.  9(b)), 
the  VLR  model  has  much  better  performance  than  other  models 
and  its  simulation  results  almost  overlap  the  experimental  data. 
The  prediction  of  the  three-branch  model  has  noticeable  discrep¬ 
ancies  from  the  measurements.  The  simulation  results  are  smaller 
than  the  measurements  during  the  initial  and  ending  phases  of  the 
simulation,  and  slightly  larger  during  the  middle  period.  The  two- 
branch  model  with  three  different  EPRs  behaves  poorly  during  the 
measured  self-discharge  period. 


2.5 


2.0 


>  1.5 

CD 

CD 

CD 

O  1.0 

> 


0.5 

0.0 

2.55 

2.50 

2.45 

>  2.40 

CD 

CD 

|  2.35 

> 

2.30 

2.25 

2.20 


a  £**AAUjJMUUUUUttUUUM 

z 

/ 

♦  experiment 
a  VLR  model 

■  three  branch 

x  two  branch  4350 

*  two  branch  2000 

•  two  branch  9000 

./ 

f 

/ 

$ 

i 

0  200  400  600  800  1000  1200  1400 

Time  (s) 


b 

. . 

. 

**  XX*xXx 

♦  experiment 
a  VLR  model 

■  three  branch 

x  two  branch  4350 

*  two  branch  2000 

*x  - ^*  %***_,  *  x  X  X  v 

**  ■■■■! 

*  X 

~**** 

*  X 

*  X 

•  two  branch  9000 

- 

0  1000  2000  3000  4000  5000  6000  7000 

Time  (s) 


Fig.  9.  Comparisons  between  different  models  of  supercapacitors:  (a)  charging  and 
redistribution  processes  and  (b)  self-discharge  process. 


5.1.  Energy  recursive  equation  (ERE)  model 

The  ERE  model  is  developed  based  on  the  leakage  power  profile 
of  the  supercapacitor.  After  a  supercapacitor  is  fully  charged,  the 
voltage  between  two  terminals  of  the  supercapacitor,  V,  is  mea¬ 
sured  during  the  leakage  process.  The  residual  energy  stored  in 
the  supercapacitor  can  be  calculated  from  the  voltage  at  different 
sampling  time  and  rated  capacitance  [5,15]. 

Residual  =  ICV2  (10) 

Then  the  energy  leaked  over  a  period  of  time  can  be  calculated. 
The  leakage  power  can  be  further  derived  by  dividing  the  leaked 
energy  by  the  sampling  period  [15]. 

The  leakage  power  profile  of  the  22F  supercapacitor  from  Cooper 
Bussmann  for  2h  is  shown  in  Fig.  10.  The  leakage  power  profile 


5.  Comparison  of  the  VLR  and  ERE  models  in  self-powered 
wireless  sensor  networks 

Supercapacitors  have  been  used  as  alternative  energy  storage 
devices  for  self-powered  wireless  sensor  networks  [3,5,12]  due  to 
their  long  charge-discharge  cycle  lives.  Accurate  supercapacitor 
models  are  critical  for  designing  efficient  energy  storage  systems 
and  supporting  power  management  related  research.  The  current 
commonly  used  models,  ERE  models,  are  based  on  leakage  char¬ 
acteristics  of  supercapacitors  [5,12].  In  this  section,  performance 
of  the  VLR  model  is  compared  with  that  of  the  ERE  model  in  a 
simplified  self-powered  sensor  node. 


Fig.  10.  Leakage  power  profile  of  a  22F  supercapacitor. 


Y.  Zhang,  H.  Yang  /  Journal  of  Power  Sources  196  ( 201 1 )  4128-4135 


4133 


can  be  approximated  with  the  following  piecewise  linear  function 
[12,15]. 


f  60.2V- 148.3,  V  =  (2.491-2.505  V) 
leakage  =<  12.61V  -27.27,  V  =  (2.434-2.491  V)  (11) 

[  6.417V- 14.63,  V  =  (2.368-2.434  V) 

Taking  the  harvested  energy  and  consumed  energy  into  account, 
the  energy  stored  in  the  supercapacitor  at  the  end  of  every  T-second 
slot  is  [12]: 


E(n  +  1)  —  E(n)  +  Eharvested(^)  ^consumed (n)  ^leakage(^)  x  T  (12) 

where  E(n  +  1)  and  E(n )  are  the  residual  energies  at  the  begin¬ 
ning  of  the  (n  + 1  )th  and  nth  time  slots,  EharvestedM  and  EconsumedM 
are  the  harvested  and  consumed  energy  during  the  nth  time  slot, 
Pieai<age(n)  is  the  leakage  power  corresponding  to  E(n)  and  is  treated 
as  a  constant  during  the  fixed  short  time  slot  T.  The  supercapacitor 
voltage  can  be  calculated  from  the  residual  energy  [5]. 


V(n  +  1)  = 


2E(n  +  l) 


(13) 


5.2.  Experiment  setup 

A  simplified  energy  storage  system  with  a  harvester  and  a  con¬ 
sumer  is  shown  in  Fig.  1 1  assuming  the  VLR  model  is  used  for  the 
supercapacitor.  If  the  ERE  model  is  used,  we  resort  to  the  energy 
recursive  equations  (Eqs.  (11 )— (13))  to  calculate  the  supercapacitor 
voltage.  The  energy  harvester  (such  as  a  solar  panel)  is  modeled  as 
a  constant  current  source  IH  [7]  with  variable  magnitude  and  duty 
cycle,  which  charges  the  supercapacitor  when  it  is  turned  on.  The 
energy  consumer  (for  example,  a  Telos  mote  [8])  is  also  modeled 
as  a  constant  current  source  IC.  It  drains  current  from  the  energy 
storage  system.  To  examine  performance  of  the  two  supercapacitor 
models,  the  energy  storage  system  is  limited  to  a  single  superca¬ 
pacitor  (the  22F  supercapacitor  used  to  illustrate  characterization 
process)  without  considering  the  input  and  output  regulators. 

The  energy  harvester  and  consumer  current  profiles  are  shown 
in  Fig.  12.  The  supercapacitor  was  completely  discharged  at  the 
beginning  of  simulation.  Then  the  harvester  (200  mA)  was  turned 
on  for  270  s,  which  charged  the  supercapacitor  close  to  its  rated 
voltage  2.5  V.  The  supercapacitor  was  subject  to  redistribution  in 


Time  (s) 

Fig.  12.  Energy  harvester  and  consumer  profiles. 
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Fig.  13.  Comparisons  between  simulation  results  of  VLR  and  ERE  models:  (a)  simu¬ 
lation  results  during  0-1700  s  and  (b)  simulation  results  during  270-1700  s. 


the  following  230  s.  From  500  to  700  s,  the  consumer  was  acti¬ 
vated  and  drained  a  constant  current  of  10  mA.  The  system  was 
idle  and  the  supercapacitor  experienced  self-discharge  in  the  next 
200  s.  Then  another  harvester  current  pulse  (10  mA)  charged  the 
supercapacitor  from  900  to  1 100  s  to  replenish  the  storage  device. 
After  that,  the  supercapacitor  repeated  redistribution,  discharging 
and  self-discharge  in  the  following  600  s. 

5.3.  Results  and  discussion 

The  supercapacitor  voltage  evolution  profiles  simulated  using 
VLR  and  ERE  models  during  the  1700  s  are  shown  in  Fig.  13(a). 
Fig.  13(b)  zooms  in  the  results  for  270-1 700  s.  In  the  first  charging 
period  (0-270  s)  in  Fig.  13(a),  the  voltages  predicted  by  two  models 
are  similar  except  at  the  end  of  this  charging  period.  The  voltage 
value  predicted  by  the  ERE  model  is  0.0959  V  more  than  that  of 
the  VLR  model  at  270  seconds.  The  voltage  difference  between  VLR 
and  ERE  models  then  increases  with  time  (Fig.  13(b)).  The  volt¬ 
age  drop  predicted  by  the  ERE  model  (0.0050  V)  is  much  smaller 
than  that  of  the  VLR  model  (0.0821  V)  during  the  first  redistribu¬ 
tion  period  (270-500  s).  This  is  because  the  VLR  model  considers 
the  internal  charge  redistribution  while  the  ERE  model  does  not.  For 
the  VLR  model,  during  the  charging  process,  not  all  internal  capac¬ 
itances  are  charged  to  the  same  voltage  level.  Part  of  the  charge 
stored  in  the  first  branch  is  redistributed  to  the  second  branch  after 
the  charging  current  is  turned  off,  which  reduces  the  voltage  value 
between  two  terminals  of  the  supercapacitor.  This  was  called  “volt¬ 
age  relaxation”  in  [32].  The  small  voltage  drop  predicted  by  the  ERE 
model  during  this  period  is  the  result  of  self-discharge,  which  is  also 
considered  in  the  VLR  model. 
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Fig.  14.  Comparison  between  experiment  and  simulation  results  of  VLR  and  ERE 
models  during  the  charging  and  redistribution  processes. 


The  accuracy  of  these  two  models  during  the  charging  and  redis¬ 
tribution  processes  can  be  evaluated  using  experimental  data  from 
Fig.  8(a).  The  simulation  setup  emulated  the  experiment.  The  con¬ 
sumer  current  source  was  turned  off.  During  the  charging  and 
redistribution  processes,  the  supercapacitor  was  first  completely 
discharged,  and  then  the  energy  harvester  was  turned  on  with  a 
constant  current  of  97  mA  for  610  s  (charging  process)  and  turned 
off  for  the  next  890  s  (redistribution  process).  The  comparison  of 
results  (Fig.  14)  shows  that  the  VLR  model  can  more  accurately  pre¬ 
dict  the  experimental  data  than  the  ERE  model  during  the  charging 
and  redistribution  processes.  The  supercapacitor  voltage  at  the  end 
of  the  charging  process  was  overestimated  by  the  ERE  model,  which 
was  also  carried  over  to  the  redistribution  process.  This  explains  the 
higher  voltage  values  predicted  by  the  ERE  model  in  Fig.  13. 

In  the  discharging  period  (500-700  s),  the  current  drained  by  the 
consumer  was  much  larger  than  redistribution  and  self-discharge 
rates.  The  voltage  drop  shown  in  Fig.  13  is  mainly  the  result  of  the 
discharge  current.  The  difference  between  the  voltage  drops  for  the 
two  models  (0.0949  V  for  ERE  model  and  0.1 241  V  for  VLR  model)  is 
still  contributed  from  the  redistribution  modeled  by  the  VLR  model. 
In  the  next  idle  period  (700-900  s),  redistribution  again  results 
in  a  larger  voltage  drop  for  the  VLR  model  (0.0408  V)  compared 
with  the  ERE  model  (0.0021  V).  For  the  second  phase  of  simu¬ 
lation  (900-1 700  s),  the  voltage  differences  between  two  models 
follow  the  same  pattern.  A  significant  voltage  difference  (0.3488  V) 
is  resulted  at  the  end  of  the  experiment  (1700  s).  The  voltage  dif¬ 
ference  between  two  models  in  Fig.  13  is  mainly  the  results  of 
the  internal  charge  redistribution.  Since  a  supercapacitor  in  a  self- 
powered  sensor  node  experiences  frequent  charging,  redistribution 
and  discharging  cycles,  the  voltage  difference  will  further  increase 
if  the  simulation  period  is  extended. 

6.  Conclusions 

An  equivalent  circuit  model  of  supercapacitors  is  pro¬ 
posed  and  the  characterization  process  is  illustrated.  A  variable 
leakage  resistance  (VLR)  is  connected  in  parallel  with  two 
resistor-capacitor  branches  to  model  the  nonlinearity  of  super¬ 
capacitor  self-discharge.  The  VLR  represents  both  leakage  due  to 
diffusion-controlled  Faradaic  redox  reactions,  which  is  a  nonlin¬ 
ear  process,  and  internal  ohmic  leakage.  The  proposed  VLR  model 
increases  model  accuracy  while  maintaining  simplicity  of  circuit 
structure.  Performance  of  the  VLR  model  is  compared  with  that  of 
three-branch  and  two-branch  models.  Results  show  that  the  VLR 
model  has  better  accuracy  in  predicting  self-discharge  of  super¬ 


capacitors  while  offering  comparable  accuracy  for  charging  and 
redistribution  processes. 

The  VLR  model  is  also  evaluated  in  a  simplified  energy  storage 
system  for  a  self-powered  wireless  sensor  node.  Performance  of  the 
VLR  model  is  compared  with  that  of  the  ERE  (energy  recursive  equa¬ 
tion)  supercapacitor  model  currently  used  for  research  in  wireless 
sensor  networks.  The  difference  between  these  two  models  is  due 
to  that  the  ERE  model  does  not  model  the  supercapacitor  internal 
charge  redistribution,  which  results  in  overestimation  of  terminal 
voltages  of  supercapacitors.  Therefore,  the  VLR  model  provides  a 
better  alternative  tool  for  supporting  research  on  energy  storage 
system  design  and  power  management  strategies  for  self-powered 
wireless  sensor  networks. 
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